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We study the structure and interfacial properties of model athermal mixtures of colloids and 
excluded volume polymers. The colloid particles are modeled as hard spheres whereas the polymer 
coils are modeled as chains formed from tangentially bonded hard spheres. Within the framework 
of the nonlocal density functional theory we study the influence of the chain length on the surface 
^L{- tension and the interfacial width. We find that the interfacial tension of the colloid-interacting 

D ' polymer mixtures increases with the chain length and is significantly smaller than that of the ideal 

polymers. For certain parameters we find oscillations on the colloid-rich parts of the density profiles 
. of both colloids and polymers with the oscillation period of the order of the colloid diameter. The 

' interfacial width is few colloid diameters wide and also increases with the chain length. We find the 

interfacial width for the end segments to be larger than that for the middle segments and this effect 
is more pronounced for longer chains. 

o 

■ I. INTRODUCTION 

' Suspensions of sterically stabilized colloids and nonadsorbing polymers often exhibit phase separation P, Q . Since 
(— I the experimental parameters can be tailored to match the desired properties, such systems have become an important 
Q tool in investigating various theoretical concepts. Recent experiments on well-characterized colloid-polymer mixtures 
O have involved studies of the wetting transition 0, 0| , real space observation of the thermal capillary waves , the 
fluid-fluid interfacial tension and the interfacial width [Tol |. 

■ One of the first theoretical models giving an insight into this phenomenon is the Asakura-Oosawa-Vrij (AOV) model 
^ . of colloid-polymer mixtures [Til IT^ , whereby the polymer-polymer interactions are ideal but the polymer-colloid and 

\^ ' colloid-colloid interactions are of the hard-sphere type. Demixing transition into the colloid-rich and the polymer-rich 
^^O phases may be explained by invoking the concept of attractive depletion interactions that arise due to a tendency 
of the system to reduce the volume excluded to the centers of polymer coils . Brader and Evans J^J have 

investigated the interfacial properties of the AOV model within the square gradient approximation. Their results were 
in a reasonable agreement with the experimental data. Subsequently Schmidt et al. |l6l | have proposed an extension 
of Rosenfeld's successful fundamental measure theory to the AOV model. It has turned out that this functional 
. yields the interfacial tension in a good agreement with the experiment a nd p redicts a sequence of layering transitions 

■ in the partial wetting regime prior to a transition to complete wetting ll8j. These findings have been confirmed by 
^ ■ Monte Carlo simulations ^3 and recent experiments 0, 3 ■ 

I ' The Asakura-Oosawa-Vrij model was originally introduced for systems, where the polymer coils are smaller than 
' ^ , the colloid dimensions and the polymer-polymer interactions are neglected. The latter approximation is reasonable 
if the second virial coefficient of the monomer-monomer interaction vanishes what corresponds to the theta-point 
conditions. The influence of the polymer-polymer interactions has been studied first by Warren et al. |2Cl| who used 
the perturbation theory around the AOV reference model. Fuchs and Schweizer |2lL |22. ,.23j have noted that the 
polymer-polymer interactions lead to an upward shift of the spinodal. These findings have been confirmed by the 
results of the modified free volume theory '23| , the geometry-based density functional approach [25!| and sophisticated 
$H computer simulation techniques .26] . Very recently the interfacial behavior of the interacting polymer-hard sphere 
5^ , colloid mixtures within the square gradient approximation has been investigated [2^ l28ll. An important ingredient of 
the square gradient approximation is an expression for the direct correlation function [29j . This is usually tackled by 
employing the mean-spherical approximation, so that only an expression for the interaction potential is needed. Both 
groups have adopted different approaches. To obtain the interaction potential between two colloidal spheres mediated 
by excluded volume polymers Aarts et al. [23| have used the generalized Gibbs adsorption equation approach [s^, 
whereas Moncho-Jorda et al. have used a depletion potential obtained from direct computer simulations |3l[. 
Both groups have found that the interfacial tension is lower for interacting polymers than for the AOV model. 
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In a recent work Paricaud et al. |33| have studied the phase behavior in model athermal mixtures of colloids 
and freely jointed tangentiall y b onded hard sphere polymers by means of Wertheim's first order thermodynamic 
perturbation theory (TPTl) They have found that if the colloid diameter is much larger than the polymer 

segment diameter the system may undergo a demixing transition. This theory offers several advantages over previously 
mentioned approaches. As a fully microscopic theory it is free of the coarse-graining, which is a powerful concept 
on its own [^^^ but may be difficult to implement, for instance, close to geometrically structured substrates [s^llsfij . 
The two-body potentials resulting from the coarse-graining techniques may be insufficient to correctly describe the 
underlying complex system and the higher-order terms may have to be taken into account 37]. Also, within the TPTl 
theory it is straightforward to incorporate the attractive interactions. However, the Wertheim TPTl description is 
not free of its own caveats. For example, an extension of the theory to the solid phase while in principle possible 
|38l |. it is difficult to carry out for the colloid-polymer mixtures. Within the TPTl theory the second virial coefficient 
scales linearly with the chain length M [S^jllfl instead of ~ M'^'^ , where v is the Flory exponent. The consequence is 
that within the TPTl treatment the dilute regime is not described accurately. Although it is possible to improve this 
deficiency it is not clear how this approach could be extended to the colloid-polymer mixtures. Despite these 
issues, the approach of Paricaud et al. 32] seems to be quite interesting. 

In the present work we apply the microscopic density functional theory of Yu and Wu 42] to study the interfacial 
properties of the hard sphere colloids and excluded volume polymers modeled as freely jointed tangentially bonded 
hard spheres. In the bulk limit this functional yields an equation of state identical with that of Ref. |32 ]. The applied 
functional is nonlocal, therefore it should correctly capture the structure of the interface including possible oscillatory 
behavior of the density profiles. Also, the Yu and Wu theory treats the polymer chains and colloid spheres on equal 
footing, therefore the many-body effects are automatically included. 



II. THEORY 



In this work we consider mixtures of colloids (species C) and polymers (species P). The colloidal particles are 
represented as hard spheres of diameter ac whereas the polymer coils are modeled as chains built of M tangentially 
bonded hard-sphere segments of diameter aps- We assume that there are no torsional or bending potentials imposed 
on the polymer segments, i.e. the monomers are freely jointed and any arbitrary polymer configuration free of the 
intermolecular and intramolecular overlaps is allowed. By introducing the total bonding potential Vb(R) as a sum of 
bonding potentials Wf, between the monomers Vf,(R) = X)f=i ^ '^b{\^i+i ~ the intramolecular interactions for the 
freely-jointed tangential hard spheres can be conveniently written as 

exp[-/3T4(R)] = 11 . (1) 

1=1 P'^ 

In the above R = (ri, r2, • • • , ym) denotes a set of monomer coordinates. 

As in every density functional theory the grand potential of the system is assumed to be a functional of the local 
densities of polymers p_p(R) and colloids jOc(r) 

f][pp(R),pc(r)] -F[pp(R),pc(r)] + 

jdRpp{K){vZ\^)~lip) 

fdrpc{r)iv}^\r) ^ ^ic) , (2) 



where vj^^CR), /ip, v}^f\r) and pc are the external and the chemical potentials for polymers and colloids, respec- 
tively. The free energy of the system F is a sum of the ideal and excess contributions, F — Fid + Fex- The ideal part 
of the free energy is known exactly 

l3F,d[pp{'R.),Pc{v)]^ P JdIlpp{R)Vb(R) 
fdRpp{K)[\n{pp{R.j) - 1] 
+ /drpc(r)[ln(Pc(r))-l], (3) 



where the dependence on the irrelevant thermal wavelengths for polymers and colloids was suppressed. 
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In their study of bulk properties of athermal mixtures of colloid and excluded- volume polymers Paricaud et al. |33| 
employed the first order thermodynamic perturbation theory of Wertheim . Within this approach the excess free 
energy of the polymer system can be treated as a sum of the excess free energy of the reference system containing un- 
bonded monomers and a perturbation due to formation of the polymer coils. Yu and Wu 42] incorporated Wertheim's 
first-order perturbation theory into the framework of the fundamental measure theory (FMT) of Rosenfeld jlTj pro- 
viding thus a convenient mean of description of inhomogeneous mixtures of tangentially jointed hard-sphere chains. 
The connection between TPTl and FMT can be made by assuming that F^x is a functional of the local density of 
colloids and average segment densities Pps{'^) defined as 

M M „ 

PPS (r) = PPS..^ ('^) ^ E / '^^^^'^ ~''^)PP (I^) ' (4) 

i—l i—1 

where /Cps,i(r) is the local density of the polymer segment i. The convolutions of the local densities Pj{r), j = PS, C 
with a suitable set of weight functions (r), a = 3, 2, 1, 0, V2, VI yield the weighted densities Ua = n^a^^ + na'^ 



n(f ) (r) = j dr'pj {r')w^^'> (r - r') . (5) 

The weight functions •* (r) are related to the geometric properties of colloidal particles and monomer segments ^3 

^^'■^(r) = e(^-|r|) , u;«(r) = <S(^-|r|), (6) 



- y^SC-f - |r|) , w[^\r) = , (7) 

4^'(r) = ^-^ , = . (8) 

In the FMT approach the excess free energy is obtained as a volume integral over the free energy density F^^ — J dr^ 
expressed as a simple function of the weighted densities. Following Yu and Wu we assume that = $p -|- ^hs, 
where ^hs describes the reference mixture of hard spheres of diameters ac and crps, while the excess free energy 
density due to the chain connectivity $p is an "inhomogeneous counterpart" of the perturbation term in TPTl. 

There are several expressions for the hard-sphere part ^hs- For the present problem we choose the elegant and 
inspiring White-Bear version of the FMT 43, 44] 

^Hs{{na}) = -noln(l - ns) H 

1-^3 

,/3 Q na (1 - 713)2 ln(l - na) 

+ {ni - 3n2nv2 ■ nv2) , .^n ^2 • ^ 

367r(n3)^(l - na)^ 

The contribution <f>p due to chain connectivity is evaluated using Wertheim's first-order perturbation theory |33ll42| 

^piM; Wr^}; aps) = ^^n^o^'\(P^^ Hynsiips; M)] , (10) 

where = 1 — riy^'^'' • Uy'^^ /{n2^^^)'^ . yns is the Boublik-Mansoori-Carnahan-Starling-Leland ^^^^ expression 

for the contact value of the radial distribution function of mixtures hard spheres 

1-7^3 4(1-713)^ 72(1 -7l3)'^ 

with ^ = 1 - ny2 • ny2/('T-2)^- 

At equilibrium the first functional derivative of the grand potential with respect to the densities of colloids and 
polymers vanishes 



J»[pp(R),Pc(r)] _ m[pp{^),pc{v)] 

6pp{n) 5pc{v) ' ^ ^ 
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which leads to the following equation for the average segment density profile 



pps(r) = exp(/3^p) / dR^S{r -Vi) exp 

i=i 

where Az(r;) is 



M 



-/314(R)-/35I^'(rO 



1=1 



(13) 



with vi{vi) being an external potential acting on the l-ih. segment. Equation can be rewritten as 

M 

/Ops(i") = exp(/3^p) ^ exp[-/3A,;(r)]G',(r)Gji/+i_j(r) , (15) 

i=l 

where the propagator function G'i(r) is determined from the recurrence relation 

G,(r) = / dv' eM-PW)]^-^^^P4—^G,_,{v') (16) 

for i = 2,3, ... ;M and with Gi(r) = 1. Likewise Eq. (|12|l leads to the following equation for the density profile of 
colloids 

pc{v) = exp [p^ic - (r) - pSF.^/Spcir)] . (17) 

Equations p5|l and (|17|l can be solved numerically using the standard Picard iterative method. 

III. RESULTS 

In homogeneous systems the thermodynamic properties can be calculated analytically. To this end we note that 
in the bulk the vector weighted densities vanish whereas the scalar weighted densities become proportional to the 
corresponding bulk densities. For the polymer weighted densities we have n^a"^^ = S,a'^Mp'p = Cc/^Vps' with 
^g'^-' = n/Gapg, ^2^-' — T^cr^g, Ci^'' = crps/2 and ^g^-* = 1. Similar, for the colloid weighted densities we have ni*^-* = 
Co^Vc'*- The total configurational free energy per unit volume, $„, can be straightforwardly obtained by inserting the 

expressions for the weighted densities into Eqs. [®-l|nil], = $H5 + $p + /3-Vp^[ln(pp^)-l]+/3"Vc^[ln(Pc')-l]- 
Finally, the pressure P, and the chemical potentials of both species are determined from 



For bulk systems our theory is identical to the Werthcim TPTl theory for athermal mixtures of hard spheres and 
excluded-volume polymers. Paricaud et al. [s^ conjectured that TPTl predicts a demixing transition into colloid-rich 
(polymer-poor) and colloid-poor (polymer-rich) phases. The coexisting equilibrium densities (binodals) were obtained 
from the condition of the equality of pressure and chemical potentials of both species in the two demixed phases. The 
spinodal lines delimiting the regions of the stability against fluctuations of density and composition were calculated 
from det[d'^^y/dpf'^dp'j^^] — 0, i,j — P, C. The critical points were evaluated from ^3 

dp^^^dp^^^ dip^V 

Figure n shows the fluid-fluid phase equilibria (solid lines), the spinodals (dashed lines) and the critical points 
(black dots) for model athermal mixtures of colloids and excluded-volume polymers considered in this work. The 
calculations were carried out for the constant colloid to polymer segment size ratio d = cc/cFps = 10 and for M = 80, 
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100, 120, 200 and 500, from top to bottom, respectively. The diagrams are plotted in the colloid packing fraction 
- reduced pressure representation, which will be useful in the presentation of the numerical results described below. 
Note that the reciprocal of the reduced pressure plays the role similar to the temperature in the case of simple fluids, 
consequently the tie lines connecting the coexisting fluid states are horizontal. The Wertheim TPTl theory leads to 
the incorrect Flory exponent {v = 0.5 rather than 0.588) therefore the comparison in terms of the gyration radius to 
colloid radius size ratio q = Rg/Rc is not entirely meaningful but on the other hand, Rg does not enter explicitly the 
numerical calculations at any point. Therefore, assuming that the polymer gyration radius scales as Rg = apM'^aps 
with ap w 0.5 41, 48, 49], our systems would correspond to q = 0.894, 1, 1.095, 1.414, and 2.236 for M = 80, 100, 
120, 200, and 500, respectively. Thus we consider both regimes q < 1 and q > 1. We note that the reduced critical 
pressure pP^'^'^^aQ and the critical colloid packing fraction rj'^^'' decrease as the chain length increases. However, as 
mentioned in Ref. js^, for M ^ cxd and for fixed d both j3P^^'^a^ and ?7c7^ tend to a finite, nonzero value. For d = 10 

we find that (3P^''^<7% = 0.725456 and r/^''^^; = 0.176477 with the leading order correction ~ 0{l/^/M). This result 
is in a qualitative agreement with the phase behavior in the so-called "protein limit" of colloid-polymer mixtures 

We turn now to the interface between the demixed fluid phases. In Figs. El and |3I we display the representative 
examples of the colloid and average polymer segment density profiles evaluated for M = 200 and for the reduced 
pressures PPa^ = 7 (A), 6 (B), 4 (C) and 3.58 (D). The profiles close to the critical point are diffuse and smooth. 
However for larger pressures oscillations develop on the colloid-rich side of the profiles (see the insets). It is interesting 
to note that the consecutive maxima of these oscillations are located at distances approximately equal the colloid 
diameter. The presence of the oscillations on both profiles can be explained by invoking the theory of the asymptotic 
decay of the correlation functions in simple fiuids [53 ■ However, here the polymer is represented as a chain consisting 
of 200 monomers and it is not clear how the arguments of Ref. '5^ can be extended to complex fluids. Nevertheless 
it is good to see that the results from the present work are qualitatively similar to the findings of Brader et al. J^] . 

Figure El shows the reduced interfacial tension j3^(jQ as a function of the order parameter Aiyc i.e. the difference 
in the colloid packing fractions in coexisting "liquid" (colloid-rich) and vapor (colloid-poor) phases evaluated for the 
chain lengths M — 80 (soHd fine), 120 (dashed fine), 200 (dotted line) and 500 (dot-dashed line). Wc find that the 
surface tension increases with the chain length. Note that the reduced surface tension is significantly smaller than 
that of ideal polymers. These findings are in accordance with the recent results obtained from the squar gradient 
approximation [27l|28|. 

In Fig. [SI we examine the critical behavior of the interfacial tension resulting from the present theory. Part a shows 
the log-log plot of the interfacial tension vs. the order parameter Aryc whereas part b displays the logarithm of the 
interfacial tension plotted against the logarithm of the relative distance to the critical point measured in terms of 
i^p _ p{cr)y p(cr) ^Yic analog of the reduced temperature r = (T^^'') - r)/T(='') for simple fiuids). We find that on 
approaching the critical point the surface tension vanishes with classical (mean field) exponents. This is not surprising 
in view of the fact that the Wertheim TPTl th eory can be derived by considering the chemical potential of a single 
chain in the reference solvent of monomers |33 . l54l |. Similar results have been obtained in a related study of binary 
polymer blends j55l |. 

The width of the interface can be quantitatively characterized by the parameter VF, defined as jH^ 

-1 

(20) 

Z = Zo 

where zq is given by pps(zq) = {\/2)[pps{z — oo)+ ps{z = — oo)] and pps(,z = oo) and pps{z = — oo) are the segment 
densities of the coexisting polymer- rich and polymer-poor phases. In Fig. [^ we plot the interfacial width calculated 
for the systems with M = 80 (solid line), 200 (dashed line) and 500 (dot-dashed line). Similar to other studies [28| . 
we find that the interfacial width increases with the chain length. For the states well removed from the critical point 
it is of the order of two colloid diameters. On approaching the critical point W diverges as expected. In simple fluids 
the interfacial width should behave similar to the correlation length ^ i.e. W ^ 5^ ^ t^^ . The inset illustrates that in 
our approach the interfacial width indeed diverges with the mean field exponent v = 1/2, what is consistent with the 
surface tension results. 

Finally we address the issue of the distribution of the particular segments across the interface. Several authors have 
pointed out that the fluid-fluid interface in polymeric systems is enriched in end-segments [5^ l58i l59j . It is related 
to a different chain length dependence of surface tension from bulk properties |63|- Figure [7| shows the middle (solid 
lines) and end segment density profiles for the systems with M — 200 and for the reduced pressures PPa^ = 7 (A), 
4.5 (E) and 3.85 (D). For the system well removed from the critical point [marked as (A)], one observes a visible 
difference between the segment profiles. The interface of the middle segments is sharper whereas the interface of the 
end segments is more diffuse. For the intermediate value of the pressure [system (E)] the difference in the distribution 
becomes less noticeable, and in the vicinity of the critical point [system (D)] the distributions of the middle and end 



W = [pps{z = oo) - pps{z = -oo) 



dppsjz) 
dz 
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segments are practically identical. This behavior is more pronounced for longer chains, as it is clearly visible in Figure 
|S1 where we plot the interfacial width of the middle and end segments evaluated for the chains with M = 80 (solid and 
dotted lines) and for the longest chains considered in this work, M = 500 (dashed and dot-dashed lines). For states 
far away from the critical point the difference of the interfacial width for M = 80 is around one polymer segment 
diameter while for M — 500 is around the colloid radius. We conclude that for still longer chains the difference in the 
distribution of the middle and end segments may be as large as colloid diameter which in turn should be significant 
since the interfacial width of the average polymer segment is of the order of 2-3 colloid diameters. This finding cannot 
be ascertained if the coarse-grained approach to the description of the polymer is employed. 

IV. CONCLUSIONS 

We have studied the interfacial properties of model athermal mixtures of colloids and excluded volume polymers. 
The colloid particles are modeled as hard spheres while the polymer coils are modeled as freely jointed tangcntially 
bonded hard sphere chains. Model systems have been studied within the framework of the nonlocal density functional 
theory for the fixed ratio of the colloid to polymer segment diameters and for different chain lengths. We have 
found that the interfacial tension of the interacting polymer-colloid mixtures increases with the chain length and is 
significantly smaller than that of the ideal polymers. These results are in accordance with the recent studies obtained 
within the square gradient approach 27, 28]. For certain parameters we find oscillations on the colloid-rich parts of 
the density profiles of both colloids and polymers with the oscillation period of the order of the colloid diameter. The 
interfacial width is few colloid diameters wide and also increases with the chain length. We find the interfacial width 
for the end segments to be larger than that for the middle segments and this effect is more pronounced for longer 
chains. Within our approach the interfacial properties close to the critical point behave in a mean field fashion. 

Before the quantitative comparison with experiments is possible some important the issues associated with the 
framework proposed by Paricaud et al. [s^ should be solved. In particular, further research should focus on tackling 
the problem of the inaccurate description of the polymer in the dilute limit and the incorporation of the fiuid-solid 
equilibria into the framework. However, despite these deficiencies the microscopic model of the colloids-excluded 
volume polymer mixtures studied in the present work seems to be an interesting alternative to the coarse-grained 
approaches. 

Acknowledgments 

This work has been supported by KBN of Poland under the Grant 1P03B03326 (years 2004-2006). 



[1] R. Tuinier, J. Rieger, and C.G. de Kruif, Adv. Colloid Interface Sci. 103, 1 (2003). 

[2] W.C.K. Poon, J. Phys.: Condens. Matter 14, R859 (2002). 

[3] W.K. Wijting, N.A.M. Besseling, and M.A. Cohen Stuart, Phys. Rev. Lett. 90, 196101 (2003). 

[4] W.K. Wijting, N.A.M. Besseling, and M.A. Cohen Stuart, J. Phys. Chem. B 107, 10565 (2003). 

[5] D.A.G.L. Aarts, M. Schmidt, and H.N.W. Lekkerkerker, Science 304, 847 (2004). 

[6] D.A.G.L. Aarts, J.H. van der Wiel, and H.N.W. Lekkerkerker, J. Phys.: Condens. Matter 15, S245 (2003). 

[7] E.H.A. de Hoog and H.N.W. Lekkerkerker, J. Phys. Chem. B 103, 5274 (1999). 

[8] B.-H. Chen, B. Payandeh, and M. Robert, Phys. Rev. E 62, 2369 (2000). 

[9] B.-H. Chen, B. Payandeh, and M. Robert, Phys. Rev. E 64, 042401 (2001). 

[10] E.H.A. de Hoog, H.N.W. Lekkerkerker, J. Schulz, and G.H. Findenegg, J. Phys. Chem. B 103, 10657 (1999). 

[11] S. Asakura and F. Oosawa, J. Chem. Phys. 22, 1255 (1954). 

[12] A. Vrij, Pure Appl. Chem. 48, 471 (1976). 

[13] A.P. Cast, C.K. Hall, and W.B. Russell, J. Colloid Interface Sci. 96, 251 (1983). 

[14] H.N.W. Lekkerkerker, W.C.K. Poon, P.N. Pusey, A. Stroobants, and PB. Warren, Europhys. Lett. 20, 559 (1992). 

[15] J.M. Brader and R. Evans, Europhys. Lett. 49, 678 (2000). 

[16] M. Schmidt, H. Lowen, J.M. Brader, and R. Evans, Phys. Rev. Lett. 85, 1934 (2000). 

[17] Y. Rosenfeld, Phys. Rev. Lett. 63, 980 (1989). 

[18] J.M. Brader, R. Evans, M. Schmidt, and H. Lowen, J. Phys.: Condens. Matter 14, LI (2002). 

[19] M. Dijkstra and R. van Roij, Phys. Rev. Lett. 89, 208303 (2002). 

[20] P.B. Warren, S.M. Ilett, and W.C.K. Poon, Phys. Rev. E 52, 5205 (1995). 

[21] M. Fuchs and K.S. Schweizer, Europhys. Lett. 51, 621 (2000). 

[22] M. Fuchs and K.S. Schweizer, Phys. Rev. E 64, 21514 (2002). 



7 



[23] M. Fuchs and K.S. Schweizer, J. Phys.: Condens. Matter 14, R239 (2002). 

[24] D.A.G.L. Aarts, R. Tuinior, and H.N.W. Lckkcrkcrkor, J. Phys.: Condens. Matter 14, 7551 (2002). 

[25] M. Schmidt, A.R. Denton, and J.M. Bradcr, J. Chcrn. Phys. 118, 1541 (2003). 

[26] P.G. Boluis, A. A. Louis, and J.-P. Hansen, Phys. Rev. Lett. 89, 128302 (2002). 

[27] D.A.G.L. Aarts, R.P.A. Dullens, H.N.W. Lekkerkerker, D. Bonn, and R. van Roij, J. Chem. Phys. 120, 1973 (2004). 

[28] A. Moncho-Jorda, B. Rotenberg, and A. A. Louis, J. Chem. Phys. 119, 12667 (2003). 

[29] R. Evans, in Fundamentals of Inhomogeneous Fluids, edited by D. Henderson (Dekker, New York, 1992), p. 85. 

[30] R. Tuinier, H.N.W. Lekkerkerker, and D.A.G.L. Aarts, Phys. Rev. E 65, 060801 (2002). 

[31] A. A. Louis, P.G. Bolhuis, E.J. Meijer, and J.-P. Hansen, J. Chem. Phys. 117, 1893 (2002). 

[32] P. Paricaud, S. Varga, and G. Jackson, J. Chem. Phys. 118, 8525 (2003). 

[33] M.S. Wertheim, J. Chem. Phys. 87, 7323 (1987). 

[34] A.A. Louis, P.G. Bolhuis, J.-P. Hansen, and E.J. Meijer, Phys. Rev. Lett. 85, 2522 (2000). 

[35] P. Bryk, R. Roth, M. Schoen, and S. Dietrich, Europhys. Lett. 63, 233 (2003). 

[36] W. Li and H.R. Ma, Eur. Phys. J. E 12, 321 (2003). 

[37] A. Moncho-Jorda, A.A. Louis, P.G. Bolhuis, and R. Roth, J. Phys.: Condens. Matter 15, S3429 (2003). 

[38] C. Vega and L.G. MacDowell, J. Chem. Phys. 114, 10411 (2001). 

[39] T. Boubhk, Mol. Phys. 68, 191 (1989). 

[40] T. Boublik, C. Vega, and M. Diaz Pena, J. Chem. Phys. 93, 730 (1990). 

[41] C. Vega, J.M. Labaig, L.G. MacDowell, and E. Sanz, J. Chem. Phys. 113, 10398 (2000). 

[42] Y.-X. Yu and J. Wu, J. Chem. Phys. 117, 2368 (2002). 

[43] R. Roth, R. Evans, A. Lang, and G. Kahl, J. Phys.: Condens. Matter 14, 12063 (2002). 

[44] Y.-X. Yu and J. Wu, J. Chem. Phys. 117, 10165 (2002). 

[45] T. Boublik, J. Chem. Phys. 53, 471 (1970). 

[46] G.A. Mansoori, N.F. Carnahan, K.E. Starling, and T.W. Leland, J. Chem. Phys. 54, 1523 (1971). 

[47] J.S. Rowlinson, Liquids and Liquid Mixtures, (Butterworths, London, 1959), p. 205. 

[48] J. Dautenhahn and C.K. Hall, Macromolecules 27, 5399 (1994). 

[49] L. Lue, J. Chem. Phys. 112, 3422 (2000). 

[50] P. Bolhuis, E.J. Meijer, and A.A. Louis, Phys. Rev. Lett. 90, 68304 (2003). 

[51] R.P. Sear, Phys. Rev. E 66, 51401 (2002). 

[52] P. Bryk, Phys. Rev. E 68, 62501 (2003). 

[53] R. Evans, R.J.F. Leote de Carvalho, J.R. Henderson, and D.C. Hoyle J. Chem. Phys. 100, 591 (1994). 

[54] A. Ben-Naim, J. Chem. Phys. 54, 1387 (1971). 

[55] P. Bryk and S. Sokolowski, J. Chem. Phys. 120, 8299 (2004). 

[56] J. Fischer and M. Mcthfcsscl, Phys. Rev. A 22, 2836 (1980). 

[57] G.J.A. Ypma, P. Cifra, E. Nies, and A.R.D. van Bergen, Macromolecules 29, 1252 (1996). 

[58] S.K. Kumar, M. Vacatello, and D.Y. Yoon, Macromolecules 23, 2189 (1990). 

[59] J. Reiter, G. Ziflferer, and O.F. Olaj, Macromolecules 23, 224 (1990). 

[60] L Szleifer and B. Widom, J. Chem. Phys. 90, 7524 (1989). 



8 




FIG. 1: The binodals (solid lines), the spinodals (dashed lines) and the critical points (black dots) for the demixing transition in 
coUoid-cxcludcd volume polymer mixtures resulting from the TPTl theory. The phase diagrams are evaluated for the constant 
ratio of colloid diameter to polymer segment diameter d = ac/crps = 10 and for M = 80, 100, 120, 200 and 500 from top to 
bottom, respectively. The dotted line marks the line of critical points of the demixing transition for different chain lengths. 




FIG. 2: Colloid density profiles across the liquid-liquid interface calculated for the system with d = 10 and M = 200. The 
profiles are evaluated for the reduced pressures /3Pal = 7 (A), 6 (B), 4 (C) and 3.58 (D). The inset shows an enlargement of 
the profile 
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FIG. 3: Average polymer segment density profiles across the liquid-liquid interface calculated for the system with d = 10 and 
M = 200. The profiles are evaluated for the reduced pressures fiPa^ = 7 (A), 6 (B), 4 (C) and 3.58 (D). 




FIG. 4: The reduced surface tension f3^ac plotted against the difference in the colloid packing fractions in the coexisting 
liquid (colloid-rich, L) and vapor (colloid-poor, V) phases, Ar]c = — Vc"'' ■ chain lengths are given in the Figure. 
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FIG. 5: Critical behavior of the reduced surface tension of coUoids with d = 10 and excluded-volurne polymers with AI — 100. 
The open circles denote the DFT results whereas the solid lines have slopes corresponding to the classical (mean-field) critical 
exponents, (a) The logarithm of the reduced surface tension /3-y(Jc plotted against the logarithm of the difference in the colloid 
packing fractions in the coexisting liquid (colloid-rich, L) and vapor (colloid-poor, V) phases, (b) The logarithm of the reduced 
surface tension P^a^ plotted as a function of the logarithm of the relative distance to the critical point in terms of the pressure 
difference, (p _ p('='-))/p(<=''). 




FIG. 6: The interfacial width as a function of the relative distance to the critical point expressed in terms of the pressure 
difference (P - P<"''')/P<"'''. The main plot shows the interfacial width calculated for M = 80 (solid line), 200 (dashed line) 
and 500 (dot-dashed line). The inset shows a double-logarithmic plot of the same quantities evaluated for M — 100. The open 
circles denote the DFT results whereas the solid lines have slopes corresponding to the classical (mean-field) critical exponents. 
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FIG. 7: Middle (the solid lines) and end (the dashed lines) segment density profiles across the liquid-liquid interface calculated 
for the system with d = 10 and M = 200. The profiles are evaluated for the reduced pressures pPa^ = 7 (A), 4.5 (E) and 3.58 
(D). 




FIG. 8: The intorfacial width of the segment density profile Wi as a function of the relative distance to the critical point 
expressed in terms of the pressure difference {P — p('^'"' )/p('^'"' . The solid and dotted lines denote the interfacial width of the 
mid (i = 40) and end [i = 80) segment density profile for 80-mer whereas the dashed and dot-dashed lines denote the interfacial 
width of the mid {i = 250) and end {i = 500) segment density profile for 500-mer. 



